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We investigate numerically charge, spin, and entanglement dynamics of two electrons confined 
in a gated semiconductor nanowire. The electrostatic coupling between electrons in the nanowire 
and the charges in the metal gates leads to a self-trapping of the electrons which results in soliton- 
like properties. We show that the interplay of an all-electrically controlled coherent transport of 
the electron solitons and of the exchange interaction can be used to realize ultrafast SWAP and 
entangling \/SWAP gates for distant spin qubits. We demonstrate that the latter gate can be 
used to generate a maximally entangled spin state of spatially separated electrons. The results 
are obtained by quantum mechanical time-dependent calculations with exact inclusion of electron- 
electron correlations. 

PACS numbers: 73.21.La, 03.67.Lx, 73.63.Nm 


Introduction. One of the most striking manifestations 
of the quantum laws of physics is entanglement m®. 
The ability to entangle qubits is also an essential ingre¬ 
dient for quantum computation [3]. The spins of elec¬ 
trons confined in an array of electrostatically defined 
quantum dots (QD) emerged as a promising candidate 
for encoding quantum bits of information SIS]. Spin 
qubits weakly interacting with their environment can 
be electrically controlled and show potential in scalabil¬ 
ity 0- Recent experiments have demonstrated fast and 
precise manipulation, initialization, and measurement of 
spin qubits confined in lateral QDs [THE] and nanowire 
QDs [151 . Furthermore long spin decoherence times have 
been reported, reaching ~ 200ns m for InAs QDs and 
~ 270/zs [15] f° r GaAs QDs. 

However, the ability to couple and entangle solid state 
spin qubits over long distances, which is essential for real¬ 
izations of scalable quantum computer architectures and 
for applications of fault-tolerant quantum error correc¬ 
tion (QEC) schemes, still seems to be one of the key 
challenges to overcome. First encouraging steps toward 
coupling remote spin qubits via microwave cavities |17l - 
[20] have been recently reported [2UH22) . Coherent long- 
range spin qubit coupling based on cotunneling phenom¬ 
ena [23] has been demonstrated for an array consisting of 
three quantum dots [24j [25]. There are also several inter¬ 
esting proposals for coupling spatially separated QD spin 
qubits, e.g. using ferromagnets [25], floating gates m , 
Majorana edge modes [251125] or superconductors [3Uj . 

Another promising platform useful for coupling spin 
qubits are mobile electrons shuttled by surface acoustic 
waves 2DE2] or flying qubits [S3], however, spin entan¬ 
glement of such moving electrons has not been reported 
so far. 

In this paper we propose an all-electrically controlled 
and ultrafast method for realization of the SWAP and 
\/SWAP gates and for generating maximally entangled 


spin states of spatially separated electrons. Our scheme 
does not require coupling with an additional external sys¬ 
tem - ‘quantum bus’ - which may simplify its implementa¬ 
tion. The proposed scheme is based on the interplay be¬ 
tween the exchange interaction and on-demand coherent 
transport of self-trapped electron solitons [34] [35] con¬ 
fined in gated semiconductor nanostructures. Exchange 
interaction, which has a short-range character, limits the 
ability to couple spatially separated spin qubits confined 
in stationary QDs. However, for mobile electron solitons 
this is not the case. Self-trapping allows for transport¬ 
ing spatially separated initially not entangled electrons 
to the region where they can entangle their spins due 
to the exchange interactionand finally be separated and 
transported back to distant regions as an entangled en¬ 
tity. To be specific we present the results for structurally 
defined InAs nanowires. However, one can expect quali¬ 
tatively similar results for electrostatically defined quan¬ 
tum wires in 2DEG/2DHG systems as proposed e.g. in 
Refs. [35], and for different materials. One can integrate 
SWAP and i/SWAP gates with all-electrically controlled 
single ultrafast quantum logic gates [37, 38; which can be 
arranged in a 2D scalable register (Ref. (39]) and be se¬ 
lectively manipulated. Such an architecture may be par¬ 
ticularly suitable for implementation of powerful QEC 
surface codes m- 


Model. We consider two electrons confined in an InAs 
semiconductor nanowire covered by 7 electrodes eL, e#, 
ej, and ei_ 4 , to which the voltages Vl, Pr, Vj, and 
Pi— 4 , are applied HU. respectively. The inter-electrode 
distance is about lOnrn. The radius of the nanowire is 
l = 5nm. The nanowire is separated from the metal by a 
dielectric material (InAlAs) 42], and the distance from 
the center of the nanowire to the metal is d = 15nm. 
The presented system can be described by the quasi-ID 
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FIG. 1. (color online) Schematic view of the nanostructure 
and its cross-section. The nanowire is covered by 7 electrodes 
which are labeled by ei _4 and €l,r,j- Two electrons with 
opposite spins (blue arrows) are confined in the nanowire. 
The positive charge induced on the gate surface is marked by 
red areas. 


Hamiltonian, 

w VVo D (. x i, X2 

(1) 

where Pj = — ihd/dxj is the x-component of the momen¬ 
tum operator for the j-th electron, j = 1,2, m = 0.023?no 
is the effective mass of the electrons confined in the InAs 
nanowire, m o is the free electron mass. The two-electron 
wave function is represented as [43] 

v{x 1 ,x 2 ,t)= ■ ( 2 ) 

It has to be antisymmetric with respect to simul¬ 
taneous exchange of the space and spin coordinates: 

= — l 4 >ji(x 2 ,xi,t), where i,j =t,4 indicate 
the spin projections of the first and the second electron 
on the z axis. 

The electrons in the nanowire induce a positive charge 
on the surface of the metal electrodes which in turn 
leads to the self-confinement of the electron wave func¬ 
tion along the wire. This wave function has soliton-like 
properties [35]: It can be transported in the form of a 
stable wave packet which maintains its shape while trav¬ 
eling. Moreover, it can reflect or pass through obstacles 
(potential barriers or wells) with 100% probability while 
preserving its shape. This effect can be exploited to real¬ 
ize on-demand transport of self-confined electrons, whose 
motions can be fully controlled by geometry and volt¬ 
ages applied to the electrodes [35]. In the general case 
(e.g. arbitrary geometry of the electrodes) the induced 
self-confining potential is determined using the Poisson- 
Schrodinger self-consistent scheme which was described 
in detail in Refs. (36l 32 [39]. However, since in the con¬ 
sidered structure electrodes form almost uniform wide 
plates (the inter-electrode distance is about lOnm), in 
order to determine the induced potential, we can use in¬ 


H — I 2^ + Vconf{xj,p(x,t )) 

.7 = 1,2 V / 


stead the image charge technique [35] . 
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This greatly simplifies and speeds up the numeri¬ 
cal calculations. Quantum calculations [44] indicate 
that it is a good approximation of the actual re¬ 
sponse potential of the electron gas. Here, e = 
14.3eo is the dielectric constant for the InAs nanowire 
[451 . The two-electron charge density is defined as 
p(x,t) = e j dx' (\^!(x,x' ,t)\ 2 + \'S’{x' ,x,t)\ 2 ), where 
\'S>(xi,X2,t )\ 2 = Xbj=t 4 \'tpij(xi,x 2 ,t)\ 2 . The voltages 
applied to the electrodes generate an additional elec¬ 
trostatic potential in the nanowire region, <pa{x, y {h z 0 ). 
which we determine by solving the Laplace equation, 
V 2 0o (x,y,z) = 0 , under the conditions <po{x,yo, z) = 
Vi , where Vi is the voltage applied to the i-th elec¬ 
trode. Thus, according to the superposition princi¬ 
ple, the total confinement potential can be expressed as 
V con f (x,p(x,t)) = -\e\V ind (x,t) - \e\(j) 0 (x,y 0 ,z 0 ). 

The electric field, E y , generated by the electrodes and 
the induced charges are also the source of the Rashba 
spin-orbit interaction in the nanowire, H R = aE y k x a z 
(46l . However, for a chosen orientation of the electrodes, 
wire, and initial electron spin (either up or down along 
the z direction) the motion of the electrons along the wire 
does not induce spin rotations. Such spin-orbit interac¬ 
tion only slightly increases spin swap times. Furthermore 
we assume that the nanowire is grown along [111] crystal¬ 
lographic direction which allows us to neglect the Dres- 
selhaus spin-orbit(DSO) interaction [45] which can affect 
slightly gate fidelity 02. 

The effective ID Coulomb interaction between charge 
carriers in a nanowire with strong parabolic confinement 
in the y and z directions has the form [45] 


(4) 

The time evolution of the system is described by the 
time-dependent Schrodinger equation ih-^'if{x\,X 2 ,t) = 
H^(xi, X 2 , t) which we solve numerically using an ex¬ 
plicit Askar-Cakmak scheme [45]. As initial condition 
we take the ground state \H(xi, X 2 , to) = ’$ l o(xi,X 2 ) for 
the two self-confined electrons under the metal elec¬ 
trodes which we obtain by solving the eigenvalue equa¬ 
tion ff'Fo(xi, X 2 ) = £^ 0 (^ 1 ) ^ 2 ) with the imaginary time 
propagation method [50]. In this approach electron- 
electron correlations are taken into account exactly. 

To characterize properties of the system we evaluate 
the spin density (i-th component) of the two-electron sys¬ 
tem, 

H f XR 

ps i {x,t) = - / dx , (4/’ l '(x, x',t)<Ji ® 

2 d XL 

+ 4 ^( 2 /, x, t)I <Ji'fy(x',x,t)), 


(5) 
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where <7,; is i-th Pauli matrix, i = x,y, z. Consequently, 
with this formula the expectation value of the electron 
spin in the left (s^) or right (s^) part of the nanostruc¬ 
ture takes the form, 

rX 0 (x R ) 

Sz {R \t)= dxp Sz (x,t), ( 6 ) 

Jx L (x 0 ) 

where Xq is the midpoint of the nanowire and xr {xr) 
is the left (right) end of the nanowire. The amount 
of entanglement between the spin in the left and the 
right part of the nanowire can be quantified by calculat¬ 
ing the concurrence m which is defined as C(p SLR ) = 
max{0, \/Ai — y/X^ — — V^t}- It varies from zero 

for completely separable (non-entangled) states to unity 
for maximally entangled states. Here, \ are eigenvalues 
(in decreasing order) of p = ps LR (& y <E> f Xy)Ps LR *( a y ® 
tj y ), ps LR is the reduced density matrix describing two- 
electron spin states in the left and right part of the 
nanowire. 

Results. First we investigate the charge dynamics and 
illustrate differences in propagation between self-trapped 
soliton-like electrons and ‘freely’ propagating electrons 
not interacting with the metal. Initially, electrons are 
confined in the nanowire under the electrodes and 
eR (see Fig. 1), which is achieved by applying Vi t 2 , 3,4 = 
— 1 meV and zero voltage to the other gates Vl,r,j = 0. 
The electrons are forced to move against each other by 
changing the voltage on electrodes e 2 and e 3 to V 2 = 
V 3=0 and by lowering the voltage on electrodes ei and 
e 4 to Vi = V 4 = —1.5 mV. The time evolution of the two- 
electron charge density p(x, t) along the wire is depicted 
in Fig. [5] It can be seen that the electrons being self- 
trapped under the metal maintain their charge density 
shape while moving, which is a characteristic feature of 
soliton waves. Furthermore, the shape is not affected by 
the collision with other electrons [Fig. [2](a)]. However this 
is not the case for ‘free’ electrons that are not interacting 
with the metal [Fig. |2jb)]. 



°0 10 20 30 40 50 60 70 80 90 100 

t(ps) 


FIG. 2. Time evolution of the charge density p(x, t ) illustrat¬ 
ing the difference in propagation and collision of two incident 
electrons between the case with (a) and without (b) image 
charge potential. Note the self-trapping and the soliton-like 
behavior in case (a). 


In order to realize a quantum \/SWAP gate using 
the nanodevice from Fig. [T] we propose the follow¬ 
ing scheme. Initially (for to) each of the electrons is 
localized under spatially separated electrodes and 
eR , respectively. Since the electrons are significantly 
away from each other (in our numerical simulation by 
about 1.2/rm) we assume that it is possible to prepare 
each of the electrons independently in a well-defined 
spin state, i.e., the electron confined under the elec¬ 
trode ei(e/j) is in spin up s^{t 0 ) = h /2 (spin down 
s^(to) = —S-/2) state. Thus in the representation (2) 
the initial state has the following form [523 : ^ 0 i x ,y) = 
( 0 ,ip L (x 1 )(pR(x 2 ),-(pR(xi)<p L (x 2 ) 1 0 ) T , where ip R (xi) 
and (pL,(%i) are the single electron ground state orbitals 
localized in the right (R) and in the left (L) dot. In 
this situation there is no entanglement between electron 
spins, i.e., C(p(to)) = 0. 


P(x,t) 

0.1 


Ps z (*,t) 
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1 

0,75 o 
0,5 (5) 
0 , 25 ^ 

0 

0 ttrap 50 100 150 200 250trel 300 

tfps) 

FIG. 3. (a) Time evolution of the charge density p(x,t) and 

scheme of the electrodes covering the nanodevice (grey). Time 
evolution of (b) the z-component of spin density ps z , (c) the 
z-component of expectation value of the electron spin in the 
left sn(t) (red) and the right srR) (blue) part of the nanowire 
(referred to the left axis). The concurrence C (black) of the 
two-electron spin state is shown with respect to the right axis 
in (c). In the inset (c’) of Fig. 3 (c) we plot the average 
distance \xi — x 2 \ between electrons for the first few collisions. 
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Then, by changing the voltage on the electrodes (in 
the same manner as in the previous example), the elec¬ 
trons start to move. When the electrons reach the re¬ 
gion under the electrode ej (for t = ttrap ~ 20ns) the 
voltage on the neighboring electrodes e 2 and e 3 is set to 
V 2 = V 3 = -lmV, which traps the electrons under ej. 
Then the electrons collide periodically under this elec¬ 
trode. The time evolution of the charge (z-th component 
of spin) density for the two electrons p(x,t ) (ps z {x,t)) is 
plotted in Fig. [3]( a) [Fig. |3|(b)]. The time evolution of the 
expectation value of the spin in the left and right part 
of the nanodevice and the concurrence C is depicted in 
Fig-K c). During each collision (due to exchange interac- 
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tion which is intrinsically present in our model) electrons 
exchange a fraction of the spin and consequently entan¬ 
glement builds up between electron spins in the left and 
right part of the nanodevice. This is illustrated in Fig. [3] 
(c’) where we also plot the average distance between elec¬ 
trons (orange line). Dips in the average value of the spin 
or concurrence is due to a local and temporary increase 
of double occupation probability [53] during the soliton 
collision. 

After many collisions for t ss 260 ps we have a situation 
where the spin density vanishes and the spin in the left 
and the right dot is equal to zero, s^(t) « sf (t) ~ 0. Fur¬ 
thermore, the concurrence reaches C(p(t)) ss 1 , which in¬ 
dicates that the spins are maximally entangled. However, 
the system is not yet in the spatially separated entangled 
state. In order to separate electrons from the area under 
the electrode ej, the voltage on the electrodes e 2 and e 3 is 
switched (for t re i ~ 260ps) to V2 = V3 = 0 , and after the 
last collision, the electrons start to move into opposite 
directions. When they reach the region under the elec¬ 
trodes and e# (initial position) for t stop ~ 280 ps, 
respectively, they are trapped again by changing the 
voltage on electrodes ei^. 3,4 to ^^, 3,4 = 2mV. Finally, 
the maximally spin-entangled state is obtained for spa¬ 
tially separated electrons characterized by the concur¬ 
rence reaching C(p) > 0.999. 

It is also instructive to analyze the probability den¬ 
sities of the components of the total two-electron wave 
function for the case when the electrons spins are not en¬ 
tangled, partially entangled, and maximally entangled, 
respectively. For the chosen initial state with opposite 
spins, during the whole evolution components with par¬ 
allel spins are zero = i/dt = 0 , while the other two 
are nonzero. The corresponding values |^M.(a;i, £2, t)\ 2 
and x 2 , t)\ 2 are depicted in Fig. H for the ini¬ 

tial moment to - non-entangled spatially separated elec¬ 
trons, t t rap - first collision, tx /2 when the concurrence be¬ 
comes C(p{ti/ 2 )) = 0.5, t re i last collision when the elec¬ 
tron spins are maximally entangled, and finally for the 
maximally entangled and separated electron spins under 
electrodes ez, and ezj. 

A similar procedure can be used to realize the two- 
qubit SWAP gate, which fully exchanges the spin of the 
two electrons. In this case electrons have to be released 
from under the electrode ej for t re i ~ 510 ps and trapped 
again under e l and en for t s t op ~ 530 ps. The results are 
shown in Fig. [ 5 ] It is clearly seen how electrons exchange 
their spins during soliton collisions under electrode ej. It 
is quite remarkable that despite many collisions the shape 
of the function (charge density) is still preserved and well 
localized. 

The proposed scheme is most sensitive to proper ad¬ 
justment of the t^ap and t s t op times. However, our 
method is quite robust against variations of t re i, the 
releasing time of electrons from under ej. For t re i — 
30 %t re i <t< t re i + 30 %t re i the concurrence reaches still 
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FIG. 4. Probability densities of the two-electron wave func- 
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panel with primes) for the instants (a) to, (b) ttrap, (c) fi /2 
[C(p(t 1 / 2 )) = 0.5], (d) t re i, and (e) t s to P during the realization 
of the \/SWAP gate. 
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FIG. 5. Same as Fig. [3] but for the SWAP gate. 


more than 90%, C(p(t)) > 0.9. In our simulations the 
voltage switching time is about 4ps long. The operation 
times r op of the proposed gates are on the order of hun¬ 
dreds of picoseconds, which is three orders of magnitude 
shorter than the reported spin decoherence time in InAs 
QDs [15]. However, one can further tune (decrease) the 
gate operation time by increasing the voltage applied to 
the gates e 3 and e 4 . 

Summary. We have shown that the interplay of 
soliton-like properties of self-trapped electrons in gated 
semiconductor nanowires and the exchange interaction 
can be exploited to realize SWAP and \/SWAI ’ gates. 
The latter gate can be used to realize maximally entan¬ 
gled spin states of spatially separated electrons. The pro¬ 
posed gates act in an ultrafast manner (subnanoseconds) 
and are controlled only by small static voltages applied 
to the electrodes which makes our proposal particularly 
suitable for addressing individual spin qubits in scalable 
quantum registers. 

We acknowledge support from the Swiss NSF, NCCR 
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